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Abstract 

55 . We derive analytical properties for the degeneracy v{N,j) occurring in the decomposition 

CN ■ ©^(-^)i)C^''^"'^ of the state space C^*^^. We also investigate the dynamics of two qubits con 

i 

pled via Ising interactions to separate spin baths, and we study the thermodynamic limit. 
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In recent years there has been an increasing interest in the description of the dynamics of 
small quantum systems interacting with their surrounding [l|. This was motivated_j3y;_ tjhe 



aaUy, 



necessity of understanding the phenomenon of decoherence in quantum systems 
and the attempt to build quantum devices that enable the implementation of quantum 
algorithms i&\. However, the main difficulty one faces in such a task consists in dealing with 
the large number of environmental degrees of freedom, which makes most of the proposed 
theoretical models impossible to be solved analytically even for finite sizes of the surrounding. 
Among the promising candidates to quantum information processing and quantum com- 
puting, spin systems seem to be the most suitable for the construction of quantum gates |7|,|8|. 
Recently, it has been shown that exact analytical solutions can be obtained for the dynamics 



ll[ |. There, the 



of few central qubits coupled to spin baths of finite and infinite sizes |9|, [lO|, 
interaction Hamiltonians together with the baths Hamiltonians are functions of the collec- 
tive spin operators of the environments. In order to derive the reduced density matrix of 
the central qubits, the partial trace over the environmental spin degrees of freedom was 
carried out within the subspaces corresponding to the different values of the total angular 
momentum of the surrounding. 

Recall that the state space of single spin-i particle is given by C^, where C denotes the 
field of complex numbers. The corresponding basis is formed by the eigenvectors {|— ), |+)} 
associated with the eigenvalues ±| of the operator S^ = ^cr^, where a^ designates the z- 
component of the Pauli operator a. In general, the state space of a system of A^ qubits is 
given by the A^-fold tensor product of the state spaces of the individual particles, namely, 
C^®^. One possible basis of the latter space consists of the state vectors 1^- |ej), with 
ej = ±. These are eigenvectors of the collective spin operator J^, where J = |X]j=i'^i- 
Alternatively, one can construct new basis composed of the common eigenvectors of the 
operators J^ and J^; we shall denote them by \j, m) such that k, < j < N/2 and —j<m< j, 
as imposed by the laws of addition of angular momentum in quantum mechanics |12|]. In 
the above, k = for A^ even, and k = 1/2 for A^ odd. Note that the scalar product of state 
vectors corresponding to different values of j vanishes. This means that the total space 
C^®^ can be decomposed as the direct sum of subspaces C^-'^^, that is 

C^'^^ = ^u{N,j)C''+\ (1) 

j = K 

The quantity v^N^j) is the multiplicity corresponding to the value j of the total angular 



moinentum; its exact form reads [1 



^^^^3)-[n/2-j) - U/2-J-J"f+J + l(f-j)!(f+j)!' ^ ^ 

Hence, given any operator G{J) on C^*^^, its trace can be written as 

N 

~2 J 

tr G = ^z/(iV,j) J2 {j,m\G\j,m). (3) 

j=K rn=-j 

Following the general ideas of the theory of open quantum systems, the problem of 
finding a relation between the multiplicities of the subspaces C^®^' and that of C^®^, where 
^^ Ni = N, naturally arises. In this work we illustrate how this problem can be solved, in 
the case A^ = A^^i + A^2, using the invariance of the trace. The latter property will also be 
used to describe the dynamics of two qubits in separate spin baths. 

A decomposition law for the degeneracy u{N,J). Let us denote by \ji,mi) the 
basis state vectors in the space C^®^' {i = 1,2). Hence the trace of G{J) can also be 
expressed as 

Ni/2 ji N2/2 J2 

trG'=^ X^ X^ X] '^(^i'ii)'^(^2,j2)(ji,j2,mi,m2|G'|ji,j2,mi,m2). (4) 

ii=Ki mi=-ji J2=K2 m2=-J2 

On the other hand we have [ij] 

jl+J2 J 

J=\ji-J2\ M=-J 

x(^' '^ Mu.M). (5) 

\mi m2 —Ml 

where the quantity in matrix form denotes Wigner 3j-symbol; obviously, the condition 
iTLi + 1712 = M along with the triangle rule |ji — J2I < J ^ ji + J2 niust be satisfied. By 
equations (jlj) and ([5]), we can write: 

junii J2,m.2 JyJ'=\h-J2\ M=-J M'=-J' 

Jl J2 J \ jl J2 J' 



v/(2J+l)(2J' + l) I \{J',M'\G\.J,M), (6) 

rrii 7712 —Ml vrrii 7712 —M' 



where we have used the fact that 3j-symbols are real. The operator G is arbitrary; it 
can be chosen such that it satisfies {J' ,M'\G\J,M) = {J, M\G\J, M)6jj>6mm' ■ In this case 
equation ([6]) reduces to 

ji,miJ2,m2 J=\jl-J2\M=-J 

(2J + 1){( ^' ^' "^ I |V,M|G|J,M). (7) 

^ \mi ra2 —Ml ^ 

The lower and upper limits of the sum over J in the above equation are, respectively, | ji — J2 1 
and ji +j2- For J < |ji — J2I, or J > ji +J2, the triple {ji,J2, J) does not satisfy the triangle 
rule and hence the corresponding Wigner 3j-symbol vanishes. Consequently, the right-hand 
side of equation ([7]) will not be affected if we take ^^'^^^ as an upper limit, and k as a lower 
limit for the sum over J such that k = for A^^i + N2 even and k = 1/2 for A^^i + N2 odd. 
This effectively allows us to exchange the order of the sums in the above equation. Then by 
comparing the resulting equation with ([3]), we obtain 

uiN, + N2,J) = J2 E i^iN,,j,MN2,j2)i-iy^''-^''^'\2J+l) 

jl,mi j2,m2 

(8) 

Herein, we have replaced M by its maximum value J (or equivalently by —J because of 
the symmetry) since the sum does not depend on this quantum number; once again the 
condition mi + 1112 = J is implied. 

Equation ([8]) can be regarded as a decomposition law for the degeneracy; many useful 
relations satisfied by the latter can be easily obtained from it. Let us first begin by noting 
that 

iV 

j^y{N,J)=(/_ Y (9) 

J=K V 2 '^/ 

JV 

^(2J+l)KiV,J)=2^. (10) 

The first equation can be readily proved by expanding the sum over J . The second one 
simply expresses the fact that the sum of the dimensions of the subspaces C^-'^^ is equal to 




the dimension of the total state space, C^*^^. Furthermore, if we let J to take the value 
^^^^^ in equation ([8]), we obtain 

jl,mi j2,m2 

Now let us suppose that J = 0, which is possible only when A^i and N2 are either both 
even or both odd positive integers. Here it should be noted that the denominator of the 
corresponding Wigner 3j-symbol contains the product (ji — J2)Kj2 — ji)' [14]; but since 
x! = 00 for X < 0, we conclude that when J = 0, the quantity under the sum sign in the 
right-hand side of equation ([8]) is nonzero only when ji = J2- In fact one should have |12l. Il4| 




By inserting the latter expression of Wigner 3j-symbol into equation (JHl), and performing 
the sum over J2 and m2, we obtain 

*- 2 ' 2 J J / -■ \2{j—m) 

j m=~j 

■ r ^1 ^21 

= ^ Z/(iVi,J>(iV2,j), (13) 

i 
where we have used the fact that Ylm=-ii^^)'^"^ ~ (~l)^"'(2j + 1)- It immediately follows 
that 

The above procedure can be easily generalized to further decompositions of the total number 
of spins. 

Dynamics of two qubits in separate spin baths. As a second application, let 
us investigate the dynamics of two qubits coupled via ising interactions to separate spin 
environments of the same size, N. The total angular momentum operators of the latter are 
denoted by J and J^. The full Hamiltonian of the composite system is given by 

H = X{alal + ay;) + ba\al + -^(^'^. + ^Uz) + /^(^l + ^') + ^Bi + ^^2- (15) 



Here, A and S are the strengths of interaction of the central qubits with each other, 7 is 
the couphng constant to the baths, and fi is the strength of an apphed magnetic field. The 
operators Hb^, with i = 1, 2, denote the Hamiltonians of the spin baths. One can show that 
the interaction Hamiltonian describing the coupling of the central qubits to the environments 
is diagonal in the standard basis of C^ (S> C^, namely. 



Hr 



7 



diag{-T,^,-A^,A^,j:, 



(16) 



where we have introduced the operators T, = J + J' and A = J — J'. Then it can be shown 
that the model Hamiltonian is given by the direct sum of the Hamiltonian operators Hi and 
H2, where 

(17a) 



7 



Hi = a,(2/i + ^S,) + UHb + 6), 
y N 



H2 = 2\a^ + -j=a,A, + UHb - 5), (17b) 

y N 

with Hb = Hb^ + Hb2 and I2 is the 2x2 unit matrix. Note that the basis vectors of the 
subspace corresponding to Hi are given by 



i) 

T) 



+ +); 



those associated with H2 are given by 



|0) =|-+), 
|i) =! + -)■ 



(18) 
(19) 

(20) 
(21) 



Thus the system under consideration can be mapped onto two pseudo two-level systems Si 
and S2 whose dynamics is governed by the operators Hi and iJ2? respectively. Each one is 
coupled to a spin environment consisted of 2N spin-^ particles with the only exception that 
Si and S2 see different compositions of the total angular momentum, namely S^ and A^, 
respectively. Notice that the above pseudo systems become completely independent from 
each other if the initial density matrix of the qubits takes the form 



P(0) 



Pll 








Pl4 





P12 


P^3 








P12 


P^3 





Ui 








pIJ 



(22) 



In such a case, it is sufficient to investigate the couphng of each pseudo system separately. 
For a reason that will become apparent bellow, we set Hb = Hb^ + Hb2 = h{Jz — Jz)-, 
where h is the strength of an applied magnetic field. Moreover, we assume that the baths 
are initially in thermal equilibrium at temperatures Ti = T2 = T (we set fc^ = 1); the 
corresponding total initial density matrix is given by 



Pb{0) = exp(-/i/3A,)/ 



2 cosh 



(23) 



-i2N 



-J- ] is the partition function. 
Under the above assumptions, the contributions of the coupling constant 6 can be neglected. 
The dynamics of S2 is quite trivial since the corresponding time evolution operator is 
diagonal. Indeed, it is easy to show that pii(t) = p^i and p44(t) = P44. Moreover, 



pi4(t) = Z-'pI Y1 E HN,j,)u{N,j,) 



Ji,mij2,m2 



(24) 



X exp{2z[2/i + 7(mi + m2)/vN]t — h(3{mi — 1112)} 
In the special case when h = or T —^ cxd, we can write 

J,M 

= pU''''oos{^y\ (25) 

For arbitrary values of h and T, the right-hand side of equation fl2^ can be evaluated within 
the computational basis; this yields 

cos2(7t/v^)-li^ 



Pl4{t)/P 



P _ 4itM 
14 — c 



1 + 



cosh2(/i/3/2) 



(26) 



Then, by expanding the cosine function in Taylor series and taking the limit N ^ 00, we 
obtain the Gaussian decay law: 

iPuit) 



P?4 



exp 



^H^ 



cosh2(/i/3/2) 



}■ 



(27) 



This means that the decoherence time scale is given by td = cosh(/i/?/2)/|7|. Obviously 
r/) — > 00 as T — i> or /i — i> 00. 

As a measure of entanglement, we use the concurrence defined by \il 

4 
C{p) = max{0, 2 maxiv^A^] - ^ v^}, (28) 

1=1 



where the quantities Aj are the eigenvalues of the operator p{ay ® <^y)p*{cry ® cry)- In our 
case, when apphed to p{t), the above definition of the concurrence leads to the evaluation 
of the eigenvalues of the operator p(t)cr^p(t)*crj. where p(t) is now restricted to the subspace 
oi Hi. A straight forward calculation yields 



C{t) = 2\puit)\. 



(29) 



An example of the evolution in time of the real value of pi4,{t) along with the concurrence 

C{t) corresponding to the initial state (| ) + I + +))/"\/2 is shown in figure [H We notice 

the revival of the concurrence in the case of finite number of spins. At short times, the 
curves corresponding to A^ ^ oo coincide with those of finite A^. 





25 



Figure 1: (Color online) Evolution in time of the real part of pii{t)/ p\^ (oscillating curve) and 

the concurrence (enveloping curve) corresponding to the initial state (| ) + I + +))/V^- Here, 

N = 100, 7 = 2, /i/3 = 1, and /i = 4. For t < 10, the curves coincide with those of the limit 

N ^ oo. 



It should be stressed that when the Hamiltonian of the composite spin bath is given by 
Hb = h{Jz + Jz) = hT,z, then 

r , , 1 27V 

Pu{t) = Pi^e'^'^' cos(7t/ViV) - i sm{jt/VN) tanh(/i/5/2) . (30) 

The existence of the sine function makes it not possible to find a relation similar to (!27|) 
when A^ ^ oo. However if we rescale the coupling constant 7 by A^ instead of viV, that 



8 



(31) 



IS ^ ^^ 

1 1_ 

exact analytical expression can be derived for the case of an infinite number of spins, namely, 



PmW =p?4exp|-zt[4/i + 7tanh(/i/?/2)]|. (32) 

Consequently the central qubits preserve their coherence, since the decoherence time scale 
in this case is infinite, as indicated by formula (l32l) . With the new scaling of 7, the larger 
the number of spins to which the qubits are coupled, the less appreciable is the decoherence. 
The Hamiltonian operator H^ can be diagonalized by dealing with the operator A^ as a 
scalar. This yields the following matrix elements in C^: 



. , X ^ sin tV4A2 + 72A2/iV 

[/22(t) = cos(tV4A2 + 72A2/Ar)+z-LA^ ^ , =^ (33) 






^ ^4A2 + ^^Al/N 

Here we have omitted the contribution of Hb = hAz since it simply introduces a global 
unitary term to the dynamics. 

Let us consider the case when the qubits are initially prepared in the maximally entangled 
state 1^) = 4j(| — h) + I H — )).( the case of the singlet state displays a similar behavior.) 
Clearly, the density matrix p(0) = \i'){ip\ belongs to the subspace corresponding to the 
Hamiltonian H2. Using the fact that |f/22(^)P + |f/23(i)P = Ib, and U22it)Ulr^{t) + U23it)Ul.^ = 
0, it can be shown that the elements of the above density matrix evolve in time according 
to p22(t) = |[1 - git)], P23 = ^[1 - fit)], where 

^ 4A7 ^ ^ A^e-^/^A.sin-(tv/4A2 + yA2/iv) ^ 

^^^ [2cosh(/i/5/2)]2^ I v^ 4A2 + 72A2/iV /' ^^ 



and 



2^2 A 2g-/./3A. sin^ (t^AX^ + 7'A^/iv) 



[2cosh(/i/3/2)]2^ I AT 4A2 + 72A2/iV 



-^^^^^A.^=====^k (37) 



AT ^4A2 + 72A2/Ar 



Figures [2] and [3] display the behavior of the concurrence as a function of time for some 
particular values of the model parameters. We can see that for hjS = 1 { i.e. at relatively 
high temperature) the concurrence shows damped oscillations and converges to a certain 
asymptotic limit which can be analytically derived, as we shall see bellow, only for h = 
and/or (3 = 0. As hf3 increases, the oscillations disappear and the concurrence converges to 
lower asymptotic values as shown in figure [2l 
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Figure 2: (Color online) Concurrence as a function of time in the case of the initial state (| — h) + 
I + -))/\/2 for N = 100, 7 = 4, /i/3 = 4, and A = 2. 



In what follows we focus our attention on the infinite temperature limit, i.e, /3 — > 0. In 
this case the reduced density matrix takes the form 

2 v - m 1 

whereas the function f{t) simplifies to 



(38) 



fit) 



2--tr^^^ 



2^2 sin^ (t^AX^ + ^^Al/N 



(39) 



N 4A2 + -f^Al/N 

Notice that < /(t) < 2 , in accordance with the general properties of density matrices in 
C^. This enables us to derive the following explicit expression for the concurrence: 



c{t) =lwm'-mt)+^ - fit)] 



(40) 
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Figure 3: (Color online) Concurrence as a function of time in the case of the initial state (| — h) + 
I + -))/\/2 for N = 100, 7 = 1, /i/3 = 1, and A = 2. 



In the thermodynamic hmit, N —>■ oo, the function f{t) can be expressed as 



m = ^Y\i- 



x^e ^^ 



■ sin^ (t^AX^ + 2-i'^xA dx. 



(41) 



4A2 + 272^2 

Some comments are in order here: We have shown in [9| that the operator Jz/vN con- 
verges to a real normal random variable a with the probability density function F{a) = 
\/2/n exp{— 2a;2}; this is also the case for the operator J'z/vN. Thus we are led to the 
task of finding the probability distribution function L{a) of the sum of two independent 
random variables ai and 02 characterized by F{ai) and ^(0^2), respectively, (note that 
the probability distribution function of aa, where a is nonzero real number, is equal to 
{l/\a\)F{a/a).) The function L{a) is simply given by the convolution of F{a) with itself, 
which yields L{a) = {l/^/^T) exp{— a^}. This becomes apparent from the change of variable 
a -^ \pl(x carried out in equation fHTl) . An other way to see that is to simply notice that 
A2/(v2iV) converges to the random variable ot ^^ F{a). From equation pTl) it follows that 

-e'Serfc(2-V (42) 

where erfc(x) denotes the complementary error function. By virtue of equation ( HOl) . we 
obtain 

C(oo) = lim C{t) = 2v^-e^$erfcf2-y (43) 

t^oo 7 V 7 / 



lim f{t) = 1 - 2v^-eVerfc(2- 
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Figure 4: (Color online) Concurrence as a function of time in the case of the initial state (| — h) + 
I + -))/\/2 for N = 100 (coincides with that of the limit iV ^ oo), 7 = 1, /i/3 = 0, and A = 2. The 
straight line corresponds to the asymptotic limit. 
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Figure 5: (Color online) C{oo) as a function of A for 7 = 2. 



In figurelHwe have plotted the concurrence as a function of time in the limit N ^ 00 along 
with the asymptotic value given by formula fHSj) . The behavior of C(oo) as a function of A 
and 7 is shown in figures [5] and [6]. As one may expect, lim C(oo) = 1, and lim C(oo) = 0. 

A^oo 7^00 
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Figure 6: (Color online) C{oo) as a function of 7 for A = 2. 



This confirms the results of 



lOj where it is shown that strong couphng between the central 



qubits reduces the effect of the environment on their dynamics. Finally it is worth mentioning 
that due to the XY interaction between the central spins, entanglement will be generated 
between them when the initial state is | ± =f). However, the corresponding off-diagonal 
elements of the reduced density matrix vanish at long times, making the asymptotic state 
of the qubits unentangled. 

In summary we have used the invariance of the trace to derive analytical properties of 
the degeneracy i/(A^, j), and to describe the dynamics of two qubits embedded in separate 
spin baths. We have shown that when the baths have the same size, the form of the model 
Hamiltonian enables us to map the full dynamics onto the evolution in time of two pseudo 
two-level systems coupled to a spin bath whose size is twice larger than the physical ones. 
This allowed us to derive the limit of an infinite number of spins within the environments 
and to analytically calculate the asymptotic state. The results of this work provide more 
evidences regarding the role played by the mutual interactions between the central qubits 
in diminishing the effects of their coupling to the surrounding spin environments. 
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